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We develop a unified model to explain the dynamics of driven one dimensional ribbon for materials with 
strain and magnetic order parameters. We show that the model equations in their most general form explain 
several results on driven magnetostrictive metallic glass ribbons such as the period doubling route to chaos 
as a function of a dc magnetic field in the presence of a sinusoidal field, the quasiperiodic route to chaos as 
a function of the sinusoidal field for a fixed dc field, and induced and suppressed chaos in the presence of an 
additional low amplitude near resonant sinusoidal field. We also investigate the influence of a low amplitude 
near resonant field on the period doubling route. The model equations also exhibit symmetry restoring crisis 
with an exponent close to unity. The model can be adopted to explain certain results on magnetoelastic beam 
and martensitic ribbon under sinusoidal driving conditions. In the latter case, we find interesting dynamics 
of a periodic one orbit switching between two equivalent wells as a function of an ac magnetic field that 
eventually makes a direct transition to chaos under resonant driving condition. The model is also applicable 
to magnetomartensites and materials with two order parameters. 
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Modeling dynamics of real physical systems is 
generally difficult. Such a task can be challeng- 
ing when a single experimental system exhibits 
a broad spectrum of dynamical features not usu- 
ally found in a single system. One such system 
was experimentally realized almost two decades 
ago in studies on driven magnetostrictive metal- 
lic glass ribbons. These investigations show that 
the system exhibits the period doubling route to 
chaos as a function of a dc magnetic field for a 
fixed sinusoidal magnetic field and the quasiperi- 
odic route to chaos as a function of the sinusoidal 
magnetic field keeping the dc field fixed. The au- 
thors reported several other dynamical features 
that includes suppressed and induced chaos in 
the presence of an additional small amplitude 
near resonant sinusoidal field, and stochastic res- 
onance. Here we design a model with strain and 
magnetization as order parameters and show that 
the model not only captures the first three fea- 
tures but also exhibits rich dynamics. The model 
equations are sufficiently general that they can 
be adopted to explain the dynamical features of 
driven magnetoelastic beam and driven marten- 
sitic ribbon. We also find interesting switching 
dynamics in the case of forced martensitic ribbon 
hitherto not reported. 
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I. INTRODUCTION 

Sinusoidally driven nonlinear oscillators have been 
used as models for a host of physical systems such as the 
charge density wave compounds 1 , Josephson junctions 2,3 , 
p — n junctions with a nonlinear element 4 , driven non- 
linear mechanical oscillators 5,6 etc. The rich chaotic 
dynamics exhibited by driven nonlinear oscillators are 
also ubiquitous to a wide variety of systems such as 
turbulence in fluid systems 7 , chemical oscillators 8 and 
cardiac tissues 9 . Since chaos is usually an undesir- 
able feature in practical applications, developments in 
chaos control 10,11 , also realized in driven nonlinear oscil- 
lators, have been applied to control the output of laser 
system 12 , enhance the performance of permanent magnet 
synchronous motor 13 and cardiac arrhythmias 14 . Natu- 
ral to multistate weakly driven oscillators in the pres- 
ence of ubiquitous noise is the stochastic resonance 15 , 
another generic noise induced cooperative phenomenon 
that enhances the signal to noise ratio. This phe- 
nomenon is also realized in several disciplines ranging 
from physics (SQUID magnetometers, optical and elec- 
tronic devises) 16 to biology 17 . 

The fact that variety of physical systems could be rep- 
resented by driven nonlinear oscillators can be ascribed 
to universality of the dynamics. However, real physical 
systems such as the magnetostrictive ribbon are much 
more complicated. While details of the material proper- 
ties may not be relevant for reasons of universality, they 
require at least two order parameters for a proper de- 
scription. Such systems abound in nature. For exam- 
ple, ferroic materials have co-existence of two or more 
order parameters such as strain, magnetization, electric 
polarization etc 18-20 . The application of a field compli- 
mentary to one order parameter (OP) induces changes in 
another OP. For example, ferroelectromagnetic materials 
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are materials that induce magnetization with the appli- 
cation of an electric field or electric polarization with 
the application of a magnetic field. Similarly, the well 
studied ferroelectric materials exhibit strain and electric 
polarization. Magnetostrictive materials and ferromag- 
netic martensites 21 also belong to the class of materials 
exhibiting strain and magnetization. Clearly, a general 
description of such systems would require at least two 
order parameters for a proper description. More than 
one order parameter also suggests that we should expect 
richer dynamics due to enhanced phase space dimension. 

Indeed rich dynamics has been realized in experiments 
on magnetostrictive metallic glass ribbon 22-27 , almost 
twenty years ago. In a series of papers, Vohra and his 
group reported a host of dynamical features in the study 
of strain bifurcations in magnetostrictive metallic glass 
ribbons subjected to the combined influence of sinusoidal 
(ac) and dc magnetic fields 22-27 . In particular, they re- 
ported the quasiperioidic (QP) route to chaos when the 
amplitude of the sinusoidal applied magnetic field h ac 
was increased in the presence of an additional dc mag- 
netic field h dc , while they found the period doubling (PD) 
route to chaos when h& c was increased keeping the ac field 
fixed 22,23 . They also reported several dynamical features 
such as the suppression and shift of the period doubling 
bifurcation point 23 , induced sub-critical bifurcation under 
small amplitude near resonant condition 24 , suppressed 
and induced chaos 25 , and stochastic resonance 27 . These 
dynamical features are observed in unannealed metallic 
glass samples unlike the PD and QP routes to chaos re- 
ported in annealed samples earlier 28 . 

Magnetostrictive metallic glasses such as METGLAS 
alloy 2605S-2(Fe78-Bi35 , ig) have striped domain struc- 
ture that lies in the plane of the ribbon. The anisotropic 
nature of magnetostrictive strain and domain movement 
under applied magnetic field leads to the softening of the 
elastic modulus that goes by the name AE effect 29-31 . 
The softening is significant even for small applied fields 
and is also nonlinear, hardening beyond a critical value 
of the magnetic field. The AE effect is also reversible. 
The buckling of the magnetostrictive metallic glass rib- 
bons under applied magnetic field (along the length of the 
ribbon) has been experimentally demonstrated 30 . Most 
of these investigations are near equilibrium studies ex- 
cept the one that exploits the nonlinear dependence of 
the Young's modulus on the magnetic field in the field- 
annealed samples 28 ' 32 where the PD and QP routes to 
chaos have been reported 28 . However, the authors em- 
phasize that rich dynamics exhibited by the metallic 
glass samples (listed above) 22-27 are not due to the AE 
effect 22,23 as the samples used are unannealed ribbons 
where the AE effect is insignificant. 

Some of these features have been explained using ap- 
propriate normal forms 33 by appealing to the universal 
nature of bifurcation phenomenon. However, since this 
approach cannot go beyond the first bifurcation, there is 
a need to develop a model that explicitly includes both 
elastic and magnetic degrees of freedom. Moreover, the 



fact that both the ac and dc fields are necessary ingredi- 
ents for most experimental features emphasizes the im- 
portance of describing magnetic degrees of freedom along 
with strain. Part of the difficulty in developing a model in 
terms of the two relevant OP's can be traced to the lack 
of knowledge about the nature of elastic and magnetic 
nonlinearities, and their mutual coupling. We recently 
developed a model set of equations for strain and mag- 
netic order parameters, and preliminary results has been 
published 34 . Here we present detailed investigations of 
the model to show that it captures a good number of 
these experimental features apart from exhibiting rich 
dynamics. 

The model equations can be adopted to describe the 
dynamics of a ribbon of a material exhibiting strain and 
magnetization under an experimental set-up similar that 
used in the study of the dynamics of a magnetoelastic 
beam 5 ' 6,35 . There are number of variants of the exper- 
imental set-up all of which have one end of the beam 
fixed and the other end is driven. (See Ref. 6 for exam- 
ple.) Here we target the results of an experiment where 
both the phase plot and the associated Poincare map has 
been reported 35 . Here a permanent magnet is attached 
to one end and the other end is driven by a low frequency 
sinusoidal magnetic field (less than a Hz) 35 . The authors 
explain the results using a parametric forcing equation 
for the strain variable. However, if the ribbon is mag- 
netoelastic, a more natural approach would be to write 
down the equations of motion for both the elastic strain 
and magnetization. As we shall see, our model equations 
explain the experimental phase plots and the Poincare 
maps. 

Our model can also be adopted to describe driven 
martensitic ribbon by omitting the magnetic degrees of 
freedom. Such experiments were performed early eighties 
in the context of internal friction studies of nonmagnetic 
martensites. The geometry is again similar to the two 
experiments considered. Experiments have been carried 
out on nonmagnetic martensite (Cu$2. 9^14. 1-^*3) sam- 
ples in the neighborhood of martensite start temperature 
T m where elastic nonlinearities are strong 36,37 . Multi- 
periodic response has been reported when the amplitude 
of the imposed vibration exceeds a critical value under 
near resonant driving conditions. In the absence of mag- 
netic degree of freedom, our equations reduce to the one 
dimensional model for martensites introduced by Bales 
and Gooding 38 . As we shall show, apart from explaining 
the multiperiodic oscillations, we find an unusual dynam- 
ics of period one orbit switching between the two equiv- 
alent wells that eventually makes a direct transition to 
chaos. 



II. A UNIFIED MODEL 

Keeping in mind that we need to develop a model that 
can be adopted to describe the three distinct physical sys- 
tems, we begin by compiling the minimal features that 
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must be included in the model. First, considering the 
geometry of the experimental set-up used in the above 
experiments, we assume a thin long ribbon (or a beam) 
of length I (larger than the other two dimensions) fixed 
at one end and set into oscillations by applying an ac 
magnetic field at the other end. Thus, for all practi- 
cal purposes, one dimensional description would be ad- 
equate. Moreover, since the displacement or strain is 
monitored at the free end, it would be adequate to use 
a single strain variable e (in one dimension) along with 
the magnetic order parameter m. Second, the elastic free 
energy is taken to have a sixth order polynomial form in 
the strain variable since such a form would also be con- 
sistent with the normal forms used for explaining certain 
dynamical features 22-26 . Moreover, a sixth order polyno- 
mial is most appropriate for describing martensites. For 
the case of magnetostrictive metallic glass samples, such 
a form can be made consistent with the speculation by 
Vohra et al 22 ~ 26 that weak magnetic nonlinearity drives 
the strong elastic nonlinearity by including appropriate 
magnetoelastic coupling. Third, as for the choice of the 
magnetoelastic coupling term, we use a weighted sum of 
symmetry preserving and symmetry breaking terms of 
the individual free energies of the two order parameters 
(for reasons explained below). Finally, keeping in mind 
that we require the general model to be reduceable to 
the case of martensites, we use a dissipative term that 
has the form of the Rayleigh dissipation function. Then, 
as we shall show later, the model reduces to the one di- 
mensional model for martensites introduced by Bales and 
Gooding 38 that supports twinning when the magnetic de- 
gree of freedom is ignored. For the case of magnetostric- 
tive ribbon, the dissipation simply serves as a damping 
mechanism. Finally, the equation for the strain variable 
is coupled to an equation for the magnetic order param- 
eter. 

The total free energy of the system has three contribu- 
tions, namely, the strain free energy, magnetic free energy 
and magnetoelastic free energy. The total free energy is 
given by 

F T = F e + F m + F me . (1) 

In one dimension, the strain variable is e = ■ Then, 
the elastic free energy is given by 

^ = Jdy'(Me(y')) + j0 2 ), (2) 

where // is the Landau free energy density and the second 
term is the gradient free energy with D being a positive 
constant that can be conveniently set to unity. Guided 
by the above considerations, we use a sixth order poly- 
nomial for the elastic free energy /; = |e 2 — §e 4 + ^e 6 . 
Here (3 and A are positive constants, but 9 can take 
on positive or negative values. The minima of the free 
energy are located at e — and e = ±e s with e s = 
[(0+ VP 2 ~ A60/A] 1 / 2 . The parameters /3, A (of the or- 
der of unity) and 9 are sufficient to control the minima. 



When 9 > 3/3 2 /4A, e = is the true minima identified 
with the high temperature phase. At 9 = 3/3 2 /4A, we 
have a first order transition with the free energy vanish- 
ing at a transformation strain e 2 = ±3/3/2 A and e = 0. 
For 9 negative, e = state is unstable. This behavior can 
be parameterized by using 9 = x f ~-^f ~ T( - wnere 2/ anc ^ 
T s are the first and second order transition temperatures 
respectively. 

The magnetic free energy is 

- m(h dc + h ac sinut)] . (3) 

Here m is the dimensionless magnetization and r c = 
T T T " with T c referring to the Curie temperature, the 
third term is the gradient free energy, hdc and h ac are 
the dc and ac magnetic fields, and to is the frequency of 
the forcing field. The constants a, b and D m are posi- 
tive. Here again, the values of the constants a and b can 
be obtained by using the magnetization at a given tem- 
perature. For all practical purposes, the three constants 
can be set to unity as has been done here. Now consider 
the magnetoelastic coupling. Generally, the nature of 
this coupling depends on the change in the crystal sym- 
metry from the parent to the transformed phase. This 
approach can not be used here as we are working in one- 
dimension. Further, in the case of metallic glass, the 
nature of the magnetoelastic free energy is not known, 
particularly since it has a glassy structure. However, 
coupling terms either preserve or break the invariance 
e — > — e and m — > —m in the elastic and magnetic free 
energies respectively. In the absence of any information, 
we model it as a weighted sum of symmetry preserving 
and breaking terms. For the metallic glass sample, the 
weighted sum is appropriate due to lack of long range 
order while the symmetry preserving term is favored for 
crystalline magnetoelastic materials. On the basis of this, 
we take the general form of the magnetoelastic coupling 
to be 

F me = -i J dy , [(l-p)e(y')m(y , )+pe 2 (y , )m 2 (y% (4) 

where £ is magnetoelastic coupling coefficient and < 
p < 1 is an adjustable weight factor. Other types of 
gradient coupling between e and m are ignored to keep 
the model simple. 

We introduce a dissipative term in the form of the 
Rayleigh dissipative function given by 

The value of the constant 7 is not easy to estimate and 
is taken to be a free parameter. While the above form 
of Fdiss merely serves as a damping term for the magne- 
tostrictive metallic glass, we emphasize that this kind of 
dissipation is essential for inducing twinned configuration 
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for martensites 38 . Indeed, the dissipative force opposing 
the abrupt motion of the twin interface is required to 
establish a mechanical equilibrium 39 . 

The kinetic energy of the system is given by 
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(6) 



where u(y,t) is the displacement variable and p is the 
density. Then, using the Lagrangian L = T — Ft and the 
Lagrange equations of motion 
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For the magnetic order parameter we use dissipative 
dynamics given by 
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Using this and the total free energy Ft, we get 
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where T is a suitable time constant. Note that the above 
equation does not represent a pure relaxational dynamics 
since m(y) is subject to ac and dc drives. Indeed, e(y) 
is being driven through m{y) that is itself subject to a 
sinusoidal field. We note that the free energies used for 
magnetization and magnetoelastic coupling are such that 
they can represent electric polarization as well. 

We shall now simplify these equations to suit the exper- 
imental set-up where one end is fixed and the other end 
is driven. We note that since the amplitude is monitored 
at the free end, the fundamental mode that is excited 
by the applied magnetic field is such that the length of 
the ribbon I is equal to half the wavelength of the funda- 
mental mode. We assume that only the dominant mode 
of vibration is supported. This is equivalent to using 
e(x,t) = A(t)sinirx/2l and m(x,t) — B(t)simrx/2l in 
Eqs.(8, 10). Transforming Eqs.(8, 10) into rescaled vari- 
ables (k = it/21,t = kt,VL = oj/k,r' = T/k) and using 
only the dominant mode we get 

r 

A(t) = - T e A(r) + 3/L4 3 (r) - ^-A 5 (r) + k 2 A{r) 

o 

+ 7 fci(r) + | (3pA(r)B 2 (r) - (1 - , (11) 



B(t) = -r' t c B{t) 



r/ ' [T) < k 2 B(r)-U(l- P )A 



2 w 2 

+ P A 2 (t)B(t)) -h ac sinnT -h dc ]. (12) 



The overdot now refers to the redefined time deriva- 
tive. Equations (11) and (12) constitute a coupled set 
of nonlinear nonautonomous ordinary differential equa- 
tions with several parameters. We shall use r e , r c , 7, T', 
and £ as free parameters. The parameters f3 and A are 
chosen appropriately for each case of the magnetostric- 
tive metallic glass ribbon, the magnetoelastic beam and 
the martensite ribbon. In addition, both h ac and h dc are 
the experimental drive parameters for the magnetostric- 
tive metallic glass while h ac is a drive parameter for the 
other two cases. Apart from this, we have no knowledge 
of the values of the parameters. 



III. DRIVEN MAGNETOSTRICTIVE METALLIC GLASS 
RIBBON 

We now use Eqs. (11, 12) to explain several dynami- 
cal features of the driven magnetostrictive metallic glass 
ribbon. In this case, the presence of both the ac and 
dc magnetic fields appears to be crucial for the observed 
features. The application of the dc field in the presence 
of the ac field has two effects. First, the vibrations are 
centered around the strain amplitude induced by the dc 
field. Noting that the strain amplitude in Eqs. (11, 12) 
corresponds to the situation where both ac and dc fields 
are present, the strain amplitude centered around the dc 
field can be described by redefining A to be the devia- 
tion from the strain amplitude A dc due to the dc field 
alone, ie., A — > A — Ad c - The value of Ad c needs to be 
supplied. Second, the free energy for the magnetic order 
parameter is tilted to one side. Further, the shape of the 
free energy is also affected due to the presence of non- 
linear magnetoelastic coupling term. In particular, we 
note that apart from the additive B term in Eq. (11), 
B 2 appears multiplicatively with A. This is equivalent 
to parametric forcing. This affects the shape of the effec- 
tive elastic free energy. Thus, we should expect (3 to be a 
function of hdc- To understand the nature of the depen- 
dence of j3 on h c i ci we first note that keeping /3 = 1 and 
hdc = 0, increasing the amplitude of the ac field leads 
to a period doubling cascade as expected from studies 
on Duffing like oscillators 40 . On the other hand, keeping 
hdc = and fixing h ac at a value where we find period 
one, and increasing (3 also leads to a period doubling se- 
quence. In contrast, keeping h ac at a value where chaos is 
seen (with (3 = 1) and increasing h^ c reverses the period 
doubling sequence reflecting the compensating relation- 
ship between /3 and hdc- This can be made quantitative 
by keeping h ac at a suitable value, say h ac = 9.2 corre- 
sponding to the onset of periodic four orbit, and increas- 
ing both hdc and /3 such that we always retain the onset 
of period four orbit. This gives a relation (3 « 1 + 10hd c - 
However, while (3 is nearly linear in hdc, the prefactor de- 
pends on the value of h ac used. For further calculations, 
we use this parameterized form of (3 in Eqs. (11, 12). For 
simplicity, we further set A — 1 for this case. 

The dynamics of the above system of equations is quite 
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FIG. 1. (a) Period doubling bifurcation as a function of hdc 
for t c = -0.2, r 6 = -1.0, fi = = 0.6,7 = 1-592, V = 
0.09, h ac = 10.5, and p = 0.32. (b) The corresponding largest 
Lyapunov exponent for the same set of parameters. 



complicated due to the presence of several parameters. 
Here we will not attempt to study the dynamics in the 
multiparameter space which is neither the focus of the 
present investigation nor desirable in such a large multi- 
parameter space. However, interesting dynamics can be 
expected when r c is such that two or three minima exist. 
We also use r c < 0. Thus, we first map-out the region of 
the parameter space where interesting dynamics is seen. 

We begin by investigating the possibility of realizing 
the period doubling route to chaos when hd c is increased 
for a fixed h ac . Here, we keep r e = —1 so that e = 
state is unstable although interesting dynamics is also 
observed for r e < 3/3 2 /4. Other parameters are fixed at 
t c = -0.2,0 = 1,7 = 1.592, T = 0.09,/ = 5,£ = 0.6 
and p = 0.32. All quantitative measures of the dynam- 
ics are obtained after discarding the first 1.5 x 10 5 time 
steps. Then, increasing hd c keeping h ac = 10.5 leads to 
the period doubling route to chaos. Bifurcation diagrams 
have been constructed by sampling the strain amplitude 
A stroboscopically. A plot of the period doubling se- 
quence is shown in Fig. 1(a). The associated largest 
Lyapunov exponent is shown in Fig. 1(b). The PD route 
is seen for a range of parameter values around the values 
used for Fig. 1. 

We now examine if our equations also exhibit the 
quasiperiodic route to chaos as a function of h ac keeping 
fixed. Here we choose the driving frequency O = 0.97, 
but similar results are obtained for a range of values of 
O less than unity. We retain most values of the param- 
eters used for the PD route, but use a higher coupling 
£ = 3.2 and a lower damping constant 7 = 0.0421. In 
this case, keeping hd c = 0.04 fixed, we sweep h ac from 2 
onwards. We find the first few harmonic frequencies (as 



found in turbulence of rotating fluid, see Ref. 7 ) in the 
region h ac — 2.5 to 3.8 beyond which quasiperiodicity is 
seen till 5.1 with the emergence of two incommensurate 
frequencies ( 1 = 1. 86252. ..,^2 = 1.79039.. and another 
O3 = Q,i — Q,2). All other dominant frequencies in the 
power spectrum can be represented by a linear combi- 
nation of the drive frequency (O = 0.97) and Oi and 
O2 . A phase plot of the torus in the (A, A, B) space 
is shown in Fig. 2(a) for h ac = 4.4. The correspond- 
ing Poincare map in the (A p ,A p ) plane is shown in Fig. 
2(b). The largest Lyapunov exponent vanishes in the re- 
gions of quasiperiodicity as can be seen from Fig. 2(d). 
Beyond h ac = 5.1 till 5.44, we find a phase locked re- 
gion arising from a small change in the frequencies par- 
ticipating in the preceding QP regime. In this region 
of periodicity, we find four dominant frequencies all of 
which are integer multiples (n = 23, 24, 25 and 28) of 
a new frequency O3 = 0.0746, which is a rational fre- 
quency close to 3 of the QP regime. (The correspond- 
ing Poincare map also confirms the region is periodic.) 
In the interval scanned from h ac — 2 — 7.5, we find three 
regions of chaos interrupted by windows of quasiperiod- 
icity and periodicity. The emergence of the first region 
of chaos can be attributed to these rational frequencies 
again becoming incommensurate. This region of chaos is 
followed by a transition to a second quasiperiodic regime 
{hdc = 5.53 — 5.64) followed by an abrupt transition to 
chaos (hdc = 5.64—5.73). The Poincare map correspond- 
ing to quasiperiodic region at h ac = 6.2, beyond the sec- 
ond region of chaos, is shown in Fig. 2(c). Interestingly, 
the dynamics prior to the first quasiperiodic region (be- 
fore the emergence of the first region of chaos) and that 
of after the second region of chaos has changed consider- 
ably as is evident from Fig. 2(b) and (c). The third QP 
and chaotic regions are separated by a periodic regime 
where we could only detect first five harmonics of the 
drive frequency since the transition to the third region of 
chaos is rather sharp. However, the power spectrum in 
the chaotic regime shows the presence several incommen- 
surate frequencies. The Poincare maps for the chaotic 
regions at h ac = 5.7 and 7.38 respectively are shown in 
Fig. 2(e) and (f). The quasiperiodic route to chaos is 
observed for a range of values of parameters around the 
values used for Fig. 2. 

We next investigate if our model predicts induced and 
suppressed chaos reported in experiments when the sys- 
tem is subjected to an additional small amplitude near 
resonant sinusoidal perturbation 25 . A number of theo- 
retical and experimental studies devoted to nonfeedback 
control of chaos use several variants of a weak harmonic 
perturbation such as parametric perturbation of suit- 
able frequency 41 , an additive forcing 3,42 and resorting 
to change in the shape of additive signal etc 43 . Some 
studies emphasize the importance of the phase factor as 
well. Fewer studies have been reported in the near reso- 
nance condition 44 as in the experiments of Vohra et aP 5 . 
However in the experiments, the authors make no refer- 
ence to the phase difference between the driving ac field 
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FIG. 2. (a) Quasiperiodic orbit in the (A, A, B) space for 
h ac = 4.4 keeping r c = -0.2, r e = -1.0,0 = 0.97,7 = 
0.0421, £ = 3.2, T' = 0.09, h dc = 0.04 and p = 0.32. (b) 
The corresponding Poincare map in the (A P ,A P ) plane. (c) 
Poincare map in the (A P ,A P ) plane for h ac — 6.2 .(d) Largest 
Lyapunov exponent as a function of h ac . (e) Poincare map 
for a chaotic orbit in the (A P ,A P ) plane for h ac — 5.7 and 
(f) Poincare map for a chaotic orbit in the (A p , A p ) plane for 
h ac = 7.38. 



and the resonant drive. For this reason we use a reso- 
nant drive of the form h(r) — h r sm£i r T, with ft r close 
to the first subharmonic il r — it/ 2 — <5 (or the second 
subharmonic) . In their experiments, the influence of the 
additional resonant field has been illustrated with respect 
to a transition from period-2 orbit to chaos as a function 
of h^c 25 - Thus, we first locate a direct transition from 
period two cycle to chaos as a function of hd c keeping 
the perturbing field h r = 0. This transition occurs at 
hdc = 0.04014. (The parameter values used are the same 
as for the period doubling route except for £ = 1.0.) A 
bifurcation diagram for the strobed strain amplitude A p 



is shown in Fig. 3(a). First we consider keeping S small. 
Keeping hd c = 0.04014, we apply the perturbing signal 
h r smfl r T with S = 10~ 6 . Even for small h r we find that 
the onset of chaos is delayed as is clear from Fig. 3(b). 
Another feature that is clear (from Fig. 3b) is that the 
amplitude of the period two orbit appears to switch be- 
tween values corresponding to the absence of resonant 
field and a slightly increased value in the presence of the 
near resonant field. The switching between the two val- 
ues is also intermittent with h r . Other changes can also 
be noted. For example, just beyond the delayed onset of 
chaos, we see a period five window. Thus, not only there 
is delayed chaos, there are some changes in the bifurca- 
tion diagram as well. The magnitude of the shift, though 
small, increases linearly with h r . 

We have also studied the influence of S on the nature of 
the dynamics as in experiments. For illustration we use 
<5 = 10~ 2 . In this case, on examining the amplitude A as 
a function of time, we find that the effect of the resonant 
sinusoidal drive manifests as a periodic modulation of 
the amplitude of the period two orbit with a period T = 
2n/S (= 2tt x 10 2 for 5 = 10~ 2 ). The amplitude of the 
modulation is proportional to h r . We find that these 
modulations persist even after long time (1 — 3 x 10 7 
time steps of At — 0.01). Thus, the associated Poincare 
map of the periodic orbit shows finite dispersion. Even 
so, a simple visual inspection of the bifurcation diagram 
displayed in Fig. 3(c) shows that the onset of chaos is 
hastened. This is confirmed by the Lypaunov exponent 
calculation. (It must be stressed here that even for small 
(5, such a periodic modulation of the amplitude of the 
periodic two orbit is seen. However, since the period of 
this modulation is large = 2ir x 10 6 and since the Poincare 
map is constructed over time duration small compared 
to the period of this modulation, the dispersion in the 
Poincare map is not visible for S = 10~ 6 .) One can also 
note that the magnitude of the shift in the onset of chaos 
(o) is significantly larger than that for the suppressed 
chaos as can be seen from Fig. 3(d). We also note that 
in our case, suppressed chaos is seen for small values of S 
while induced chaos is seen for large values of <5, opposite 
of what is reported. 

Identifying fio with the onset of chaos (hdc = 0.04014) 
when h r = and \i that for finite h r , the normalized shift 
(□) for the onset of chaos (obtained by computing 

the Lyapunov exponent) as function of h r is shown in 
Fig. 3(d). Using fl r = ft/2 + 5 does not affect induced or 
delayed chaos, but the magnitude of the shift is different. 

We now examine how the near resonant forcing affects 
the period doubling sequence as a function hdc shown 
in Fig. la. (The parameter values used are the same 
as for the PD route shown in Fig. 1). Consider reso- 
nant forcing with respect to the first subharmonic, ie., 
Vt r = ft/2 — 6 keeping S = 10 -6 . In the absence of h r , 
while the PD route starts with period one cycle (see Fig. 
la), even for very small resonant forcing amplitude, say 
h r = 0.02, we do not find the period one cycle. A plot 
of the bifurcation diagram keeping h r = 0.02 is shown in 
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FIG. 3. (a) Bifurcation diagram as a function of hdc for r c = 
-0.2, r e = -l.O,^ = 1-0, A = 1,7 = 1-592, T' = 0.09, h ac = 
10.5 and p = 0.32. (b) Bifurcation diagram as a function of 
hdc for the same set of parameters for S = 10 -6 and h r = 0.05 
that clearly shows delayed chaos. Note that the period two 
orbit appears to switch by a small amount intermittently, (c) 
Bifurcation diagram as a function of hdc for the same set of 
parameters for S — 10~ 2 and h r = 0.05 that clearly shows 
induced chaos, (d) Suppressed and induced chaos for 8 = 
10" 6 (D) and 5 = 10~ 2 (o) as a function of h r . 



Fig. 4(a). As can be seen from the figure, not only the 
bifurcation digram starts with period two cycle but also 
the entire bifurcation diagram is affected (compare Fig. 
1(a)). While successive PD bifurcations leading to chaos 
are hastened, a period four orbit interrupts the chaotic 
band. As we increase h^c, the bifurcation sequence also 
culminates in chaos. For larger values of ft. r (say 0.035), 
the first chaotic band is interrupted by a period two orbit 
instead of period four orbit seen for smaller h r . This is 
shown in Fig. 4(b). This feature of altering the entire 
bifurcation diagram in the presence of a resonant per- 
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FIG. 4. (a) Bifurcation diagram as a function of hdc for r c = 
-0.2, T e = -1.0, £ = 0.6, n = 1,7 = 1.592, F' = 0.09, h ac = 
10.5, p — 0.32 and h r = 0.02. (b) Bifurcation diagram as a 
function of hdc for r c = — 0.2, r e = — 1.0, £ = 0.6, f2 = 1,7 = 
1.592, F' = 0.09, h ac = 10.5, p = 0.32 and h r = 0.035. 



turbing field is similar to earlier studies on the effect of 
nonresonant perturbing sinusoidal field 44 . 

However, using larger 5, say 5 = 10~ 3 , even for small 
h r = 0.02, transients persist even after very long runs 
(10 7 time steps). The influence on period two cycle is so 
significant that one sees a band in the Poincare map in- 
stead of two points. This is due to fact that the sinusoidal 
modulations of period 2n/S with amplitude proportional 
to h r persist even after long runs. The magnitude of fluc- 
tuations become significant near bifurcation points. We 
have verified that a similar effect is also seen in the case 
of the Duffing's Oscillator. 

As is well known, driven oscillators exhibit several in- 
teresting dynamical features such as persistent chaotic 
transients noted near crisis points 45 ' 46 . For example, for 
a system with two wells, symmetry restoring crisis is usu- 
ally seen 46,47 . Thus, symmetry restoring crisis is a dis- 
tinct possibility since our model has left-right symmetry. 
Here we investigate this possibility in the case of the pe- 
riod doubling route to chaos observed as a function of 
h ac keeping hd c fixed at a suitable value. For a given 
set of initial conditions located in one well, the attractor 
is confined to the same well for a range of values of h ac 
starting from zero. However, as we increase h ac further 
the attractor suddenly expands to occupy the symmetry 
allowed phase space corresponding to the two wells at 
a critical value h ac = h c ac = 11.194. A plot of the bi- 
furcation diagram showing the sudden expansion of the 
attractor is displayed in Fig. 5(a). Note also that near 
h c ac , the attractor expands slowly to occupy larger posi- 
tive values of the original well before jumping to the left 
well. Very long runs were required near the transition 
regime. To find transient chaos as we approach the criti- 
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FIG. 5. (a) Symmetry restoring crisis during period doubling 
bifurcation as a function of h ac - The parameter values are 
the same as used in Fig. 1, except that we keep hdc = 0.04 
and /3 = 1.4 and vary h ac . Other parameter values are r c = 
-0.2, r e = -1.0, fi = 1,£ = 0.6,7 = 1-592, T' = 0.09 and 
p — 0.32. (b) Power law dependence of the mean residence 
time in the pre-attractor. 

cal value from above, we have numerically calculated the 
mean residence time < r > in the pre-crisis attractor for 
various values of h ac . We find that < r > scales as a 
power of h ac - h c ac , ie., < r (h ac - h^ c )^ n where r\ 
is an exponent 46 ' 47 . A log- log plot of the mean residence 
time in the pre-crisis attractor is shown in Fig. 5(b). The 
value of the exponent obtained turns out to be 77 ~ 1.06. 



A. Driven magnetoelastic beam 

The above model equations can be easily adopted to 
describe the experimental results of the dynamics of a 
magnetoelastic beam 5 ' 6 ' 35 . Here, we target experiments 
where both the phase plots and the associated Poincare 
maps have been reported 35 . These authors first obtain 
the static bifurcations of the beam exhibiting two equiv- 
alent symmetric strains either through subcritical or su- 
percritical bifurcation depending on the distance of the 
drive from the epoxied magnet. Consider the case of su- 
percritical static bifurcation which implies r e < in our 
model. Further, the material is a crystalline alloy which 
corresponds to p = 1 in our model. The frequency of the 
drive used in their study is much less than the natural 
frequency of the beam. Further the damping is stated to 
be low which implies that we should use small 7 value in 
our model. The authors monitor the phase plot in the 
(A, A) plane and also obtain the corresponding Poincare 
map. Our purpose is to see if these two results can be 



FIG. 6. (a) Phase plot in the (A, A) plane for r c = —0.1, T e = 
-0.25, fi = 0.13, £ = 3.7,7 = 0.05, V = 0.1, p = 1.0, h dc = 
0,/3 = A = 1 and ho = 0.42. (b) The corresponding Poincare 
plot. 



reproduced by using our model equations. 

We have numerically solved Eqs. (11, 12) keeping the 
parameters fixed at t c = —0.1, r e = —0.25, £1 = 0.13, £ = 
3.7,7 = 0.05, T' = 0.1, hdc = 0,/? = 1, and A = 1. 
Since experiments use low frequency, for illustration we 
use n = 0.13, which is nearly a factor ten smaller than 
the natural frequency of the system(£l = 1). Initially for 
small h ac we find that the strain amplitude decays. As 
we increase h ac beyond a critical value, we see a (A, A) 
phase plot similar to Fig. 5(a) of Ref. 35 . However, in 
our case, since the phase plot is a projection of the three 
dimensional flow (corresponding to elastic and magnetic 
degrees of freedom) on the (A, A) plane, the Poincare 
map can be sensitive to the value of h ac . Varying h ac , 
we find a Poincare map similar to Fig. 5(b) of Ref. 35 
for h ac — 0.42. The phase plot in (A, A) plane and the 
associated Poincare map are shown in Fig. 6(a) and (b) 
respectively. Clearly these plots are very similar to those 
in Ref. 35 . Similar Poincare maps (as in experiments) are 
seen for a range of values beyond h ac = 0.38. 



IV. DRIVEN MARTENSITES 

As stated in the introduction, the model equations can 
also be adopted to describe the dynamics of nonmagnetic 
martensitic (Cm§2. 9^14. 1^3) ribbon carried out in the 
context of internal friction experiments quite early 36,37 . 
These experiments have been performed by fixing the 
sample at one end and applying a sinusoidal drive at the 
other end. Multiperiodic oscillations have been reported. 
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Before proceeding further, we note that it is traditional 
to use a sixth order Landau polynomial in e such that the 
martensite minima occur at e = ±1. This is equivalent 
to using (3 = 2 and A = 3 in the free energy. Further, we 
note that the sample has a low Ni content, which how- 
ever, is adequate to enforce periodic oscillations using an 
applied magnetic field h ac . Alternately, small permanent 
magnet is attached that help to drive the free end of the 
sample 37 . Since the alloy is nonmagnetic, we can drop 
the equation for m and terms containing B in Eq. (11). 
Then we have an equation for the strain variable that 
reduces to the form derived by Bales and Gooding 38 for 
a martensite in one dimension. We note that the equa- 
tion for the strain variable with the dissipative term gives 
twinned solutions (near square wave type of solutions, see 
Fig. 1 of Ref. 38 ). Thus, we have a single equation for the 
strain amplitude with a driving term h ac sinVLt given by 



Mr) 



Here, r„ 



15 



(r« + k 2 )A + 6A 3 (t) - —A 5 (t) + jkA(r) 



iQt 



(13) 



T-T s 



i T -, where T m is the martensite onset tem- 
perature. Note that the effective transformation temper- 
ature is r c + k 2 . 

First we address the presence of multiperioidic oscilla- 
tions when the read is driven at near resonant frequency 
keeping the temperature just above the martensite trans- 
formation temperature T m . The resonance curve corre- 
sponding to Eq. (13) is given by 



((T e + k 2 )-n 2 )A + 



9A 3 
~2~ 



75A 5 



64 



^mA)=hl c . 

A plot of the resonance curve is shown in Fig. 7 for 
driving amplitude h ac = 1.46. Clearly there are regions 
of multistability that are realized when the initial con- 
ditions are in suitable basins of attraction. The gap be- 
tween the upper two lines closes for small values of h ac . 
Also shown in the plot are the numerically obtained res- 
onance curve marked by open circles that exhibit several 
resonances for ft < 1. 

First, we briefly summarize the results of the internal 
friction experiments 36,37 . A cursory look at the figure 
of relaxational oscillations (Fig. 2 of Ref. 36 ) shows that 
there are four distinct peaks with a rough periodicity 
of 310 s. (The nature of the curve for low amplitudes 
is considerably rough which can be attributed to lack 
of control in measuring devises or instrumental errors.) 
The authors state that these oscillations are seen beyond 
a certain amplitude of the drive. While the authors rec- 
ognize that these are nonlinear oscillations, they fail to 
recognize that they are subharmonics since they specifi- 
cally state that they are harmonics. (However, it is not 
even clear if these are stabilized oscillations.) Similarly, 
the plots of the strain amplitude given in Figure 3 and 
4 of Ref. 37 also appear to be multiperiodic. In the ab- 
sence of reliable plots of these relaxational oscillations we 
regard them as multiperiodic. 
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FIG. 7. (a) Resonance curve given by Eq. (14) for ho = 
1.46,7 — 0.25 and r e = 0.95. Also shown is the numerically 
obtained resonance curve that shows several resonances for Q 
less than unity, (b) A typical amplitude plot for = 0.99995. 



We now examine the response of Eq. (13) under con- 
ditions close to those in experiments. These experiments 
have been conducted at a temperature close to marten- 
site transformation temperature (1° above T m ) and under 
resonant forcing conditions. The natural frequency of the 
ribbon in our model is unity. We use a drive frequency 
fl = 0.99995 (that roughly corresponds to 5uj — 0.005s -1 
when the natural frequency is ~ 100). Using r e = 0.95, 
the effective temperature r e + k 2 = 1.05 corresponding 
to a temperature to just above T m . When the ampli- 
tude of the drive h ac is small, we see a wave form that 
follows that of the drive. This is related to the small 
gap between the two upper branches which opens-up as 
we increase h ac . This happens around h ac = 1.35. For 
values beyond this we see a period two cycle. A typical 
period two cycle for h ac ~ 1.46 and = 0.99995 is shown 
in Fig. 7(b). As can be seen, the wave form is qualita- 
tively similar to the experimental plot even though they 
are of period two. (See Fig. 2 of Ref. 36 .) This also corre- 
sponds to the resonance near f2 ~ 1. We also find period 
three, five and seven cycles at the resonances found with 
decreasing £1 However, we do not see period four. 

The dynamics is more interesting when the temper- 
ature is lower than Tm- Usually resonant driving con- 
ditions are expected to exhibit only periodic response. 
However, in this case, we find unusual dynamics. To il- 
lustrate this, we choose r e = — 1 with fl = 0.99995. For 
small values of h ac , with initial conditions in one of the 
wells, the strain amplitude remains confined to the well. 
Beyond a certain value of h ac , typically the barrier height 
of the potential well (~ 0.035 for these parameter values), 
the system still settles to a periodic state confined to the 
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FIG. 8. (a) Bifurcation diagram for r e = — l,fi = 
0.99995, P = 2; A = 3; 7 = 0.25 (b) A typical Poincare map 
in the chaotic regime for h ac = 0.64. 



same well. As we increase the drive amplitude h ac fur- 
ther, the transients become long. It is beyond this value 
that we find that the orbit can settle into either of the two 
wells starting from the same initial condition in one of the 
wells. The switching dynamics is seen when h ac exceeds 
0.038. A bifurcation diagram of the strobed amplitude 
A p is shown in Fig. 8a. The inset clearly shows the dy- 
namics of the period one cycle switching between the two 
wells. We could not find any regularity in this switching 
process as a function h ac (see the inset). As can be seen 
from Fig. 8a, we find period five and period one are inter- 
spersed in the region 0.057 < h ac < 058. Intriguingly, we 
could not find any pattern in the duration of these tran- 
sients as we increased h ac . Sometimes even a very small 
change in h ac led to very short transients and at other 
times it resulted in long transients. Another interesting 
feature is that for certain values of h ac , we find long tran- 
sients that appeared to be period five, which eventually 
settled to period one cycle confined to one of the wells. 
The switching dynamics is seen till h ac — 0.646. Further 
increase in h ac leads to a direct transition from the pe- 
riodic one orbit to chaos as shown in Fig. 8. We note 
here that a direct transition from a period one cycle to 
chaos has been realized in the study of Josephson junc- 
tions under nonresonant drive conditions 48 . Moreover, 
the amplitude of the periodic one orbit does not fluctu- 
ate. In contrast, in our case, the dynamics is observed 
under near resonant drive conditions. To the best of our 
knowledge, we are not aware of any reports on this kind 
of switching phenomenon. More work is needed to under- 
stand the underlying mechanism leading to the switching 
dynamics. 



V. DISCUSSION 

In summary, we have derived a set of model equa- 
tions that are general enough to be applicable to sev- 
eral driven systems with strain and magnetization as the 
order parameters. For instance, they capture a host of 
experimental results on the dynamics of driven magne- 
tostrictive metallic glass ribbon 22-26 . Two other physical 
systems considered are the dynamics of driven magne- 
toelastic beam and martensitic ribbon. In all the cases 
considered, the samples are thin long ribbons fixed at 
one end and driven at the other end. Since in each case, 
the geometry of the experimental set-up is essentially the 
same, a natural starting point is to derive a general set 
of equations that can be reduced to each of these cases. 
This is what has been achieved here. We have demon- 
strated that these set of equations do in fact recover a 
large number of experimental results. 

In their most general form, these equations explain sev- 
eral experimental results on the strain dynamics of the 
magnetostrictive metallic glass ribbons reported almost 
two decades ago 22-26 . Detailed numerical studies of the 
model equations show that these equations reproduce the 
period doubling route to chaos as a function of hd c keep- 
ing h ac fixed and also the quasiperiodic route to chaos as 
a function of h ac for a fixed dc field. The model equations 
predict induced and suppressed chaos when the system 
is subjected to an additional small amplitude near reso- 
nant perturbation with a frequency fl r close to the first 
subharmonic. In our case, the suppressed chaos is seen 
for small 8 [8 = 0/2 — O r ) while induced chaos is seen for 
large 8, opposite of what is reported in experiments. This 
may be attributed to parametric type of forcing in our 
equations. These equations also exhibit symmetry restor- 
ing crisis with an exponent close to unity. We have also 
examined the influence of near resonant forcing on the 
period doubling sequence as a function hd c - We find that 
the entire bifurcation diagram is affected even for small 
8 consistent with the earlier report 44 . While the succes- 
sive PD bifurcations leading to chaos are hastened, the 
nature of bifurcations for larger values of hd c are signif- 
icantly altered in detail. However, for large 8, the effect 
of resonant forcing is to induce modulation of the am- 
plitude of the original signal. For example the periodic 
orbits are modulated with a period 2tt/8 inducing large 
dispersion in the Poincare maps. 

The model equations have been adopted to explain the 
phase plot and the associated Poincare map of a magnc- 
toelastic beam driven at a low frequency 35 . We note 
here that these authors represent the coupling between 
magnetization and strain variables by providing a para- 
metric drive for the strain variable, instead of describing 
the magnetic degrees of freedom along with the elastic 
degrees of freedom as has been done here. Clearly, our 
set of equations are better suited to explain the exper- 
imental results as our model includes both strain and 
magnetization. 

The model when reduced to describe the driven 
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martensitic ribbon also explains the multiperiodic non- 
linear response reported in the internal friction experi- 
ments on samples of nonmagnetic martensite samples of 
Cus2.9Al14.1Ni3 36 ' 37 . It is pertinent to note here that 
our equations automatically describe the twinned struc- 
ture in one dimension when the magnetic degree of free- 
dom is dropped 38 . More interestingly, the equation for 
the strain amplitude exhibits some unusual switching dy- 
namics between two equivalent wells of the periodic one 
orbit, which eventually makes a direct transition from 
period one cycle to chaos. 

Clearly, these equations can be easily adopted to de- 
scribe the dynamics of driven magnetomartensites as 
well 21 . Indeed, rich dynamics is predicted in magne- 
tomartensites where both elastic and magnetic nonlin- 
earities are significant if experiments are performed in a 
similar geometry. While no such experiment has been 
reported so far, it would be interesting to verify this pre- 
diction. We hope that the present work will encourage 
such dynamical studies. 

The approach is clearly applicable to any class of ma- 
terials with two order parameters 18 . For example these 
equations should describe the dynamics of driven ferro- 
electric materials that can be described by strain and 
electric polarization with the latter replacing magneti- 
zation. Such systems are also described by polynomial 
forms of free energies as used in the model. Our analysis 
suggests rich dynamics if experiments are performed in a 
similar geometry on material with two order parameters. 

Here it is worthwhile to point out that Eq. (8) can 
be adopted to describe an order parameter where iner- 
tia is unimportant by taking 7 large. Then, it is clear 
that our approach can be easily extended to the case of 
ferroelectromagentic materials (such as RMnO^ where 
R = Y,In) ls . We hope that our work would encour- 
age dynamic experiments in such materials where exper- 
iments are traditionally carried out in quasistatic condi- 
tions. 
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